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Abstract 

A variational method is studied based on the minimum of energy variance. 
The method is tested on exactly soluble problems in quantum mechanics, 
and is shown to be a useful tool whenever the properties of states are more 
relevant than the eigenvalues. In quantum field theory the method provides 
a consistent second order extension of the gaussian effective potential. 
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I. INTRODUCTION 



In 1873 Lord Rayleigh [1] described a variational method for calculating the frequencies 
of mechanical systems. Since then the Rayleigh-Ritz method has been an important tool 
for the approximate solution of physical problems. In quantum mechanics the method has 
proven very useful, and provides an upper bound for the ground state energy, as the exact 
eigenstate of the hamiltonian H yields the lowest energy expectation value. 

More generally, the exact eigcnstatcs arc known to be stationary points for the expecta- 
tion value of H. That is not a special property of H: for any real function /, the expectation 
value of the operator f{H) can be proven to be stationary at the exact eigenstates of H. 
Thus, when the search of approximate eigenstates of H is the main issue, the expectation 
value of any function f{H) can be used as a functional of the trial eigenstate. In general 
the result is different and depends on the choice of the function / unless the trial state is 
the exact eigenstate of H. This dependence is a measure of the accuracy of the approximate 
eigenstates, and can be used as a variational method of calculation whenever the description 
of states is more important than the determination of the corresponding energies. 

In this paper the variance of H is shown to be the natural choice for a measure of 
the dependence on /. The outcoming variational method, which we call "Minimal Energy 
Variance" (MEV), is not novel, having been used since 1955 in numerical calculations [2]. 
The method has not been very popular as the average of the square H"^ is required, and only 
recently some of its interesting properties have been shown in numerical Quantum Monte 
Carlo calculations [3,4]. 

We discuss in some detail the properties of MEV and show that it can be regarded as 
a useful complementary tool for the properties of the eigensates more than a substitute for 
the usual variational method. In that respect the MEV turns out to be of interest even for 
analytical calculations despite the larger amount of work required for its evaluation which is 
comparable to a second order perturbative approximation. We show that MEV is at least as 
general as the standard variational method, and that it can be relavant for the variational 
treatment of quantum field theories like the scalar theory. 

We begin by defining the MEV and discussing its general properties in Sec.II. In fact 
MEV is usually described as a numerical tool in the framework of Quantum Monte Carlo 
calculations, despite its generality. In Sec. Ill the method is illustrated by comparison of some 
results for exactly solvable problems in quantum mechanics (the harmonic oscillator and the 
hydrogen atom). In Sec. IV MEV is shown to be relevant for the variational treatment of a 
scalar theory in order to get a consistent second order extension of the Gaussian Effective 
Potential. 



II. DEFINITION AND PROPERTIES 



Denoting by |^) a generic state in the Hilbert space, the expectation value of f{H) reads 



and the stationary condition is 
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which is satisfied if 

= (3) 

The solutions of the eigenvalue problem are the stationary points of (/). Usually (with 
some special exception [5]) these solutions are the eigenstates of H. Thus whatever is /, 
the stationary points of (/) yield approximate eigenstates of H. However the stationary 
point does depend on the choice of the function / if the generic trial state |\E') belongs 
to a sub-space which does not contain the exact eigenstate. For instance it is well known 
that if the trial state |^) is not an exact eigenstate then in general (i?") 7^ {H)"'. Thus 
there is no reason why the simple choice f{H) = H should be the best choice. Actually 
any different choice for the function / would give a different weighting of the trial state 
in the expectation value. While the usual energy variation seems to be reasonable for the 
approximate evaluation of the ground state energy, that might not be the best choice for 
describing the properties of the eigenstates. 

Wc assume that the trial state |^) is closer to the exact eigenstate when its sensitivity 
to / is lower. In other words we define a distance D between the trial state and the exact 
eigenstate according to 

D={fiH))-fi{H)). (4) 

Thus any choice of the function / would provide a different variational method as we know 
that D = for any exact eigenstate. The most simple non trivial choice for / is f{H) = H"^ 
and yields the variance cr^ = {H"^) — {HY as a viable candidate for the distance D. 
A formal proof that the exact eigenstates of H are stationary points of cr^ is trivial: 

(^1^)^ = [^'1^) - {H^m] - m) [Hm - mn (5) 

and the right hand side vanishes if |^) is an eigenvector of H. Moreover o"^ > 0, and it 
vanishes for any exact eigenvector so that has a minimum at any eigenvector (not just 
the ground state). In practice, whenever a trial state is close enough to an eigenstate, 
the variance is expected to show a local minimum. The value of at the minimum 
is a measure of the accuracy of the corresponding approximate eigenstate. Moreover the 
minimum of cr^ acquires a deeper physical meaning if related to the dynamical properties 
of the state. The vanishing of a for eigenstates can be seen as a consequence of the time 
delocalization of the stationary states. According to Heisenberg relations At f» /i/cr, a 
smaller energy variance allows for a longer survival of the approximate eigenstate. Thus 
MEV yields approximate eigenstates that best resemble the exact ones in their dynamical 
evolution. In that respect MEV seems to be a complementary tool for the properties of the 
cigcnsates more than a substitute for the usual variational method which always gives the 
best approximation for the eigenvalues. 

III. ANALYTICAL TESTS 

The method can be tested on exactly solvable problems: the hydrogen atom and the 
harmonic oscillator. 
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A. The hydrogen atom 



In atomic units the hamiltonian of the hydrogen atom is 

H^-\V^-1 (6) 

We choose a two-parameter trial state 

{r\^>) = N{l~ar)e-^'' (7) 

This is the exact ground state for « = 0, = 1 while it is the first excited state for 
a = P = 0.5. The expectation values of H and are easily evaluated 

{H) = [ri(a, (3) - V,{a, (5)] / [2D{a, (5)] (8) 



,„2. T,{a,P) + V,{a,P)-S{a,P) 
^""^^ ^Di^) 

where Ti = (3\a^ - a/J + /S^), Vi = l3{3a'^ - Aa/3 + 2/3^), T2 = (3\3a^ + bap + 5p^), 
V2 = 4/52 (a2 _ 2af3 + 2f3^), S = 2(3^ {a'^ + 6f3^) and D = {3a^ - 3af3 + f3^). The variance cr^ 
has two local minima for (a, P) equal to (0, 1) and (0.5, 0.5), while the energy (H) only has 
a saddle point at (0.5,0.5). For a = P the two methods yield the same result: a minimum 
for a = P = 0.5 where the trial state becomes the first excited state. For a = —P, the trial 
state is quite bad: its behaviour for r — > is 

(r|*) ~ A^(l - -p^r^ + 0{r^)) ~ A^e"^^''^' (10) 

The variance is quite sensitive to the shape of the wave function, and in this case it 
fails to show any minimum, while the energy (H) still has a minimum for P = 1.5. It is 
instructive to study the behaviour of for a constant ratio a/P = k. The trial state l^f) 
may get very close to an exact eigenstate if A; or A; 1. In Fig.l the variance cr^ is shown 
for several values of k. We observe a cross-over from a pronounced minimum at P — 0.5 (for 
A; = 1) to a pronounced minimum at P — 1 (for A; = 0). As A; gets negative and moves away 
from 0, the minimum value of (j^ raises and eventually the minimum disappears as the trial 
state becomes worse and worse. Conversely as k approaches 0, a minimum around P = I 
deepens until vanishes for A; = 0. The minumum at /3 = is always present as in that 
limit the trial wave function becomes a constant which is an exact unbounded eigenstate 
with a vanishing energy. Thus at variance with the standard variational method, MEV may 
be used for approximating excited states without having to insert orthogonality conditions: 
a local minimum appears whenever the trial state is close enough to an exact eigenstate. 
Moreover the sensitivity of discards bad approximations as the minimum disappears for 
the worse trial states. Whenever a minimum is present its value is by itself a measure of the 
accuracy of the state as cr^ = for the exact eigenstate. 
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B. The harmonic oscillator 



Other insights on the method come from the study of the simple harmonic oscillator. 
Let us consider the Hamiltonian 

^ = y + (11) 

which describes an oscillator whose frequency is ujq. Let us denote by |n, o^o) the exact eigen- 
states with energies En — hujQ{n + 1/2). As a trial state we may take a linear combination 
of the lower energy eigenstates of a generic oscillator whose frequency is uj: 

\^) = \{),uj) + a\l,uj) (12) 

where both a and uo are variational parameters. The trial state is the exact ground state 
of i7 for a = and ui = ujq, while it gives the first excited state for a — > cxd and u) = ujq. 
The calculation of the expectation values of H and if^ is trivial: in units of TujJq/2 we may 
express them as 

(i/) = cosh(lnx)/3(a) (13) 



{H^) = 3sinh2(lnx)/5(a) + /9(a) (14) 

wheie X = ujq/u) and /„(«) = (l + nQ;^)/(l + Q;^) is a smooth increasing function of a ranging 
from 1 (at q; = 0) to n (for a — > 00). 
The variance follows as 

— T^l + ^(a)sinh2(lna;) (15) 

cosh (mo;) 

where g{a) = 2(3a'' + 6a^ + 1)/(1 + a^)^ is a smooth increasing function of a ranging from 
2 (at q; = 0) to 6 (for a — > 00). 

First of all we mention that both methods must predict the exact values of a even for 
X ^ 1 (i.e. UJ ^ 

cjq)' i'^ fact at any u the states \0,uj) and \l,uj) have different symmetry 
properties, and thus the trial state can be an eigenstate of parity only for a = (even) or 
a — > 00 (odd). Actually we may observe that a — and a — > 00 are stationary points 
of {H) and cr^ for any choice of the parameter x. This is evident for (H) as in Eq.(13) 
the contributions of a and x are in different factors. We always get a minimum for a — 
(ground state), while the limit a — cxd is a maximum (first excited state). Whatever is a, 
(H) has a unique stationary point for x = 1 where the hyperbolic cosine has a minimum. 
From Eq.(15) we see that the variance follows the same path: for any choice of a a minimum 
occurs at X — 1 where the hyperbolic sine vanishes. If we set x 7^ 1, then we can explore 
the dependence of the variance on a. As we move from x = 1 (i.e. ou — ujq), the trial state 
gets worse and worse. We still find two minima at a = and a ^ 00, but the minimum 
value of the variance increases as the state gets worse. In Fig. 2 the variance is shown for 
some values of x. From Eq.(15) we see that at the minima = g(a) sinh^(lna;) so that we 
get a larger variance cr^ = 6 sinh^ (In x) at the first excited state {a — > 00), and a smaller 
variance cr^ = 2sinh^(lnx) at the ground state {a = 0). Thus the trial state is a better 
approximation for the ground state than it is for the first excited state. 
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IV. SCALAR THEORY 



In quantum field theory the properties of vacuum are more relevant than its energy 
which is not finite anyway. For instance the symmetry breaking mechanism and the mass 
of the Higgs boson depend on the structure of the true vacuum. Thus we argue that the 
use of MEV could give rise to new insights on the ground state properties of relevant field 
theories like the scalar theory. In fact we show that MEV may be used for improving the 
Gaussian Effective Potential (GEP), a useful variational tool which has been discussed by 
several authors since 1974 [6-13]. The GEP has many merits, and has been successfully 
applied to physical problems ranging from electroweak symmetry breaking [14] and scalar 
theories [13], to superconductivity in bulk materials [15] and films [16]. A second order 
extension of the gaussian approximation would be desirable for a better understanding of 
the symmetry breaking transition. In fact sometimes the GEP is known to predict a first 
order transition even when the phase change should be continuous. Moreover, the GEP 
fails to show a minimum for some ranges of parameters. Attempts to improve the GEP 
have not been so successful: the Post Gaussian Effective Potential (PGEP) discussed by 
Stancu and Stevenson [17] fails to reach a minimum for any finite value of the variational 
parameter which is fixed by the vanishing of the second derivative [18]. A way out has been 
studied by Tedesco and Cea [19] who take the variational parameter fixed at the first order 
value. In this paper we point out that the minimum of variance would be a viable tool for 
determining the variational parameter, and we show that this choice allows a useful second 
order extension of the GEP. 

The GEP can be seen as an improved first order perturbative approximation. Let us 
decompose the hamiltonian in two parts as if = if^ + where Hq, is any solvable hamil- 
tonian which depends on the parameter Q, while Vq — H — Hq. The decomposition itself 
depends on the parameter Q. The ground state of Hq satisfies the eigenvalue equation 



Then the first order perturbative approximation for the lower eigenvalue of H follows 



The minimum of E can be found by a variation of the parameter Vl, and at the minimum 
point fl — Qq we get the best decomposition of H (in the sense that the first order per- 
turbative approximation yields the lower energy). However E is the expectation value of 
the full hamiltonian H, and the method is a genuine variational method: the trial state is 
the eigenstate |\1'q) which depends on the parameter Q according to Eq.(16). In the GEP 
Hq is the hamiltonian of a free scalar field whose mass is fl, and its ground state is a 
gaussian functional of fields. 

The PGEP [17] is equivalent to the second order perturbative evaluation of the vacuum 
ground state energy (effective potential). It arises from the sum of all the second order 
connected one-particle irreducible diagrams without external legs. It can be proven to be 
equivalent to the cumulant expansion discussed by Klcincrt [20] , and then the second order 
correction 5E^'^^ is basically equivalent to the variance up to a sign 



(16) 



E^E^ + (*n|Vh|*n) 



(17) 



(18) 
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Thus the minimum of the variance is equivalent to the minimum absolute value of the second 
order correction. According to the asymptotic convergence of the perturbative expansion 
we know that a minimum of the second order correction is equivalent to a minimum of the 
error that we expect in the first order expansion. Prom this point of view the minimum of 
the variance singles out the best perturbative expansion. 

The explicit expression for the second order effective potential V^^'' has been reported 
in Ref. [17] as a function of the vacuum expectation value of the field (0) = for a scalar 
theory whose action reads 



S[(j)] = J (fx -0(x) (-9^ + m^) (t){x) + A0^(a;) 
in a d-dimensional Euclidean space. The second order effective potential is given by 



y(2)^y(i)+5£;(2) 

where V^^^ is the first order GEP 

y(^) = h{n) + JmV^ + Ay'^ + \lo{^) W - + 12A(^2 + 6A/o(l^) 
and the second order correction reads 



(19) 



(20) 



(21) 



5^(2) = _ |i/(2)(0) p - ^2 ^ 12AV92 + 12A/o(l^)]' + SAV'/^'H^^) + ^A2/W(0)| . (22) 



Here /("^(Q) and /n(^) are the integrals defined according to [17] 



(23) 



I d'^x[G{x)r (25) 
where G{x) is the free particle Green function 

Most of these integrals are diverging and must be regularized. The search for the min- 
imum of yields a gap equation for the free field mass Vt. A numerical analysis of this 
gap equation shows that there is no minimum for the second order effective potential, while 
the second order correction by itself (the variance) has a pronounced minimum for a broad 
range of the parameters. 

The second order correction has been evaluated as a function of the bare parameters m 
and A, and the variational parameter O. An energy cut-off A has been inserted in order to 
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regularize all the diverging integrals. In Fig.3 SE^"^^ versus Q is reported, for = 3 and for 
the set of parameters m^/A^ = —0.06, A/A = 0.1, (f/^/X = 0.1. According to Eq.(18) it 
turns out to be negative; moreover, as it is clear from the figure, the second order correction, 
while owning a maximum when its absolute value is minimum (which is the minimum of the 
variance: see the Fig.3 inset ), is not bounded: this explains why the total effective potential 
fails to reach a minimum for any choice of the free mass fl. 

The minimum of o"^ yields a best value 1] = for each value of the shift (f. Insertion 
in V^"^^ gives our second order effective potential. This should be compared to the PGEP of 
Ref. [17] where the best fl is obtained by the vanishing of the second derivative of V^"^^ [18]. 

In Fig. 4 our second order effective potential is reported (the same d — 3 and bare pa- 
rameter values as those in Fig.3 were used). For this set the system is close to its transition 
point. For comparison in the same figure we also show the standard first order GEP, and the 
PGEP evaluated according to Ref. [17]. For d = 3 the system may be regarded as a static 
statistical model for a phase transition in the three-dimensional space (Ginzburg-Landau ac- 
tion) . The predictions of this model can be tested by comparison with the experimental data 
on the phase transition of different systems like superfiuids and superconductors. Unfortu- 
nately the simple GEP predicts a first order transition in this case (while the transition is 
known to be continuous). In Fig. 5 an enlargement of the if = area makes these reasonings 
more evident: the GEP (dotted line) is an increasing function up to a maximum (the point 
= is a local minimum). Actually the phase transition occurs when the true minimum 
rises more than the local minimum (first order transition). Our second order potential (solid 
line), evaluated by MEV, predicts a continuous transition (as it should be), with the point 
at (f = always being a local maximum in the broken phase. Thus the method provides a 
consistent second order extension of the GEP while retaining its variational character. 

We conclude that while MEV has been recently shown to be a useful tool in numeri- 
cal Quantum Monte Carlo calculations, its potentialities have not been fully explored yet. 
Whenever the properties of states are more relevant than the eigenvalues, MEV provides a 
viable variational method which can be used in analytical and field theory calculations as a 
complementary tool. 
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FIGURES 

FIG. 1. The variance for approximate eigenstates of the hydrogen atom. The trial wave 
function is defined according to Eq.(7) with a = kf3 and A; = 1.0, 0.8, 0.6, 0.4, 0.2, 0.0, -0.2, -0.4, 
-0.6, -0.8. The minimum moves from left to right when k decreases, and disappears at /c ~ —0.7 

FIG. 2. The variance o"^ for approximate eigenstates of the harmonic oscillator. The trial state 
is defined according to Eq.(12) with a ranging from a = (ground state) to a — oo (first excited 
state) while uq/oj = x is taken to be x = 1 (lower curve), x = 1.25 and x = 1.5 (upper curve). 
Approaching the exact eigenstates {x 1) the minima decrease and eventually vanish. 

FIG. 3. The second order correction (the variance, up to a sign) for the effective potential 

in scalar theory. The values for the parameters are: cZ = 3, m?' / h? = —0.06 and A/A = 0.1, 

tp/s/X = 0.1. Note that it is unbounded,so that the total effective potential cannot have a minimum 

whatever the variational parameter Q is; as shown in the inset graph (an enlargement with 

unit scaled by a factor 10^), however, it has, by itself, a pronounced maximum, making MEV a 

reasonable alternative to PGEP, where the vanishing of V^'^^ second derivative is required. 

FIG. 4. The second order effective potential evaluated by the method of minimum variance 

(sohd line) for d = 3, rr? j h? = -0.06 and A/A = 0.1. For comparison the PGEP (dashed line) 
and the simple first order GEP (dotted line) are reported. The effective potential is scaled by a 
factor 10^, while the field shift ip is in units of \/A. 

FIG. 5. An enlargement (the effective potential is scaled by a factor 10^) of the <^ = region 
in Fig.4; the GEP (dotted line) predicts a first order transition with the point at <^ = being a 
local minimum; instead, the MEV (solid line) predicts a continuous transition (as it should be for 
the superconductivity), thus providing a consistent variational second order extension of the GEP. 
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